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ABSTRACT 

We present estimates for the accretion rates in 13 X-ray bursting sources which 
exhibit photospheric expansion, basing on theoretical models of stationary, 
radiatively driven winds from neutron stars. The relatively high values 
obtained, M acc J> 10~ 9 M@/yr, are in accordance with theoretical limits for 
unstable helium burning, and, at the same time, almost never exceed the 
"dynamical" limit for stationary accretion, ~ IOMem- The only exceptions are 
1820-30, already known to be a very peculiar object, and 1608-522; there are 
indications, however, that in both sources, accretion could be non-stationary. 
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I. INTRODUCTION 

Low-Mass X-ray Binaries (LMXBs) are widely believed to be powered by the 
accretion of gas lost by the secondary, main sequence component of the binary 
system onto a neutron star. This appears to be now a well established point 
on both observational and theoretical grounds, although the observed X-ray 
properties of LMXBs can differ substantially from source to source, according 
to the different physical conditions in the system. The presence of a high 
neutron star magnetic field can produce a regular modulation in the X-ray 
emission, as in X-ray pulsators, while other subclasses of sources, in the large 
family of LMXBs, may be characterized by different accretion scenarios. Within 
this picture, one of the most important parameters in governing the overall 
appearance of the source, is certainly the mass transfer rate, M acc , from the 
secondary onto the neutron star. It has been recently suggested (Kylafis & 
Xilouris 1993) that Super Soft Sources ultraviolet emission is produced in the 
very dense, optically thick envelope which shrouds a neutron star accreting 
close to the maximum allowed rate. In this respect, X-ray bursters, too, 
should be characterized by quite high values of the accretion rate since unstable 
nuclear burning, which is responsible for the bursting phenomenon, can take 
place only if M acc exceeds some definite limit (Fushiki & Lamb 1987; Taam 
et al. 1993). The determination of accretion rates in X-ray bursting sources 
with photospheric expansion, based on independent arguments, would be quite 
useful both for testing the consistency of the helium burning scenario, and for 
shedding light on various aspects of evolution of close binaries. Any direct 
measurement of mass transfer rates in LMXBs appears to be still beyond the 
capabilities of present instrumentation and, up to now, the only, indirect, 
estimates come from the comparison of the binary period variations with 
existing evolutionary theories (see, for instance, Rappaport et al. 1987). In the 
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present paper we derive accretion rates in LMXBs which show strong X-ray 
bursts with photospheric expansion. During the expansion/contraction phase, 
in fact, a supersonic wind is believed to be present and the envelope mass at 
the beginning of the wind phase can be obtained by confronting the observed 
spectral parameters with a set of theoretical wind models (Nobili, Turolla & 
Lapidus 1993); the knowledge of the interburst time will then provide the 
accretion rate. The analysis of the 13 sources for which sufficient data are 
available, produced M acc > 10 -9 M & /yr, in agreement with the current idea 
that strong X-ray bursters should be characterized by high values of the mass 
transfer rate in comparison with other LMXBs. 

II. PHOTOSPHERIC EXPANSION IN X-RAY BURSTERS 

The most powerful observed X-ray bursts show a characteristic temporal 
behavior. Usually the event starts with a sudden increase of the X-ray intensity, 
with rise times less than one second, followed by a decrease, with a total duration 
of a few seconds. After the precursor, a noticeable decay of the flux is observed, 
lasting up to ~ 10 s in the strongest bursts, at the end of which the luminosity 
can be as low as the persistent one. The main part of the burst then begins. 
The increase of the X-ray intensity first appears in the soft energy channels, 
and gradually becomes visible in the harder bands. The color temperature 
increases while the X-ray luminosity stays nearly constant at its maximum 
value, commonly associated with the Eddington limit. A direct use of the 
relation 4TraR^ ol T^ ol = L ~ LEdd shows that the typical size of the emitting 
region, R co i, after having reached a maximum, decreases in this phase until R co \ 
is, again, comparable with the neutron stars radius R*. After the blackbody 
temperature has reached a maximum (often above ~ 2 keV), the decay starts 
with the progressive decrease of the X-ray flux accompanied by a softening of 
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the spectrum, while the "color" radius remains approximately constant. The 
last phase is quite similar to that one observed in weaker, type I, X-ray bursts, 
where the energy released is below the Eddington limit during all the burst and 
no expansion is observed. 

Photospheric expansion during strong bursts is widely believed to be 
produced by a supersonic outflow, driven by the large, super-Eddington 
luminosity released by thermonuclear He burning at the base of the envelope. 
In a very recent paper Nobili, Turolla & Lapidus (1993, hereafter paper I, see 
also Lewin, van Paradijs & Taam 1993, hereafter LVPT, for a review of earlier 
studies on this subject) have presented a more complete model for radiative wind 
acceleration during strong, type I, X-ray bursts which accounts properly for 
both energy production by 3-a helium burning and Comptonization heating- 
cooling in the outflowing envelope. One of the major results obtained in that 
investigation was the discovery of a lower limit for the mass loss rate M m i n , 
below which no stationary, supersonic wind can exist. The existence of such a 
bound is due to a sort of "preheating" effect and could be of great significance 
for the determination of some important physical parameters in strong X-ray 
bursting sources. In paper I, in fact, we suggested that the quasi-stationary 
expansion/contraction phase, during which L ~ LEdd, could be thought as a 
sequence of steady wind models with decreasing M which terminates when 
M m i n is reached. By taking this into account and if the maximum and 
minimum color temperatures in the expansion phase, T™f x , T™l n , are known 
from observational data, the initial envelope mass M env can be derived (see 
paper I for details); the mass accretion rate M acc follows once the time interval 
At between two successive bursts with photospheric expansion is available. The 
fitting of the theoretical R co i -T co \ relation with the the observed one provides 
also an estimate of the spectral hardening factor 7. In paper I this approach has 
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been applied to 4U/MXB 1820-30, mainly as a working example of the method. 

Photospheric expansion during powerful bursts has been observed so far in 
14 sources of type I bursts, and in type II bursts from the Rapid Burster 1730- 
335 (LVPT), although it can not be excluded that in some cases the physical 
scenario may differ from that one outlined above. Here we present results for 
all the bursters for which published data allow a derivation of T^ x , T™ ; m and 
At. 

III. INDIVIDUAL SOURCES 

In this section we shall review the main observational properties of all sources of 
type I bursts in which photospheric expansion has been detected. Our primary 
goal is to extract from observational data the maximum and minimum values of 
the color temperature during the bursts with radius expansion together with the 
interburst time. As we discussed in the previous section, in fact, if these three 
parameters are known both the hardening factor and the accretion rate onto 
the bursting neutron star can be derived from the sequence of stationary wind 
models. The set of models we shall use below refers to a neutron star of mass 
M* = 1.5M and radius i?* = 13.5 km. For 1820-30 there is a strong evidence 
that the secondary is an evolved, helium-rich star, so helium wind models are 
used, while, in comparing data relative to all other sources, "solar" composition 
models seem more appropriate. For each source we derive 7, M acc and compare 
the adiabatic cooling time of the envelope, t coo i, with the interburst time. As 
it was discussed in paper I, accretion can be regarded as stationary only if 

tcool <§C At. 

a) 0748-676 

During 9 observations with EXOSAT a total of 37 bursts were detected 
(Gottwald et al. 1986). Bursts with expansion always had At > 5 hr, and they 
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occurred only when the persistent flux was high, ~2x 10~ 9 ergcm~ 2 s _1 . In 
the paper by Gottwald et al. the data of all bursts are mixed together, and 
only average figures of interest for us may be extracted: T™ ; ax « 2.2 keV, and 

T min ^ 1A ke y Then 7 = 1-6g and Menv = 4 x 1Q 22 g? with aR accretion rate Q f 

M acc ~ 3.5 x 10~ 8 (5hr/ At)M Q /yr. The usual assumption that the persistent 
flux is due to the conversion of gravitational potential energy appears to be 
justified because the cooling time of the envelope is t coo \ ~4x 10 3 s < At and, 
therefore, accretion is stationary, contrary to what occurs in 1820-30 (see paper 
I and section IV). 

Gottwald et al. (1987) observed also the source in a low state when 
no expansion occurred during bursts. These observations seem to confirm 
the suggestion by Fujimoto et al. (1987) that incomplete nuclear burning is 
responsible for bursts with short At, while in bursts with expansion (which 
show larger At) nearly all the fuel is burned out. 

b) 1 516-569/ Cir X-l 

EXOSAT observations (Tennant, Fabian, & Shafer 1986) revealed just one 
burst with photospheric expansion. The spectral data give T™]* x m 2.2 keV, 
T™\ n R3 1.5 keV, which implies 7 = 1.68, M env fh 3 x 10 22 g. The time interval 
separating this burst from the previous one is unknown, but it can be estimated 
~ 4 days. The estimated accretion rate M acc & 1.4 x 10~ 9 M & /yr is below the 
critical value and accretion should be stationary since t coo \ ~ 2.5 x 10 3 s <C At. 

c) 1608-522 

The observational data of Hakucho (Murakami et al. 1987) contain the 
record of a burst with radial expansion that occurred on April 8, 1980. The 
bursts frequency in April-May was i/& = 0.09 ± 0.03 hr" 1 . The blackbody fitting 
of spectra gives T™ ax ss 3.5 - 4.0 keV, T™ m w 1.2 keV. The hardening factor 
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is 7 = 2.67, and a quite large value of the envelope mass is obtained, M env = 
2 x 10 23 g. The high cooling time of such an envelope, t coo \ ~ 10 5 s J> At, and 
the high, largely super-Eddington accretion rate M acc ~8x 10~ 8 M Q /yr make 
the source, at that period of time, more similar to 1820-30 than to the "normal" 
bursters with M acc ~ Msdd- 

Later observations by Tenma in 1984 indicate a decrease of the accretion 
rate in this system. Nakamura et al. (1989) reported a series of bursts, two of 
which (G and J in their notation) showed radial expansion. In burst G the fitting 
temperatures were T™? x = 2.5 keV, T™\ n = 1.8 keV, resulting in 7 = 1.91 and 
M env = 2.5 x 10 22 g. In the period of time preceeding burst G, the average 
accretion rate was M acc -6x 10~ 9 M Q /yr. Burst J was stronger, with color 
temperatures T c ™ ax = 2.9 keV and T™j n = 1.5 keV implying 7 = 2.21, M env = 
8 x 10 22 g and a super-Eddington accretion rate M acc ~ 1.7 x 10~ 8 M Q /yr. 
The accretion regime was, probably, stationary in both bursts since t coo \ ~ 
1.7 x 10 3 s < At = 6.7 x 10 4 s (burst G), and t cool ~ 2 x 10 4 s < At = 7.4 x 10 4 
s (burst J). 

Photospheric second-range oscillations, observed in 1608-522 during the 
long (~ 12 s) flat top of the light curve, are discussed in a separate paper 
(Lapidus, Nobili & Turolla 1994). 

d) 1636-536 

This object is one of most "reliable" sources of bursts and has been 
extensively observed by Hakucho, Tenma and EXOSAT. As far as Tenma 
observations are concerned, there is only one burst, denoted as burst H by 
Inoue et al. (1984), which data can enable us to determine M env and M acc . 
In all the other 11 bursts they reported, either there is no expansion or (for 
bursts D and E) the value At before the burst with expansion is unknown since 
observations were not continuous. Burst H had T™? x ~ 2.9 keV, T™{ n ~ 2 keV 
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which correspond to M env ~ 3.2 x 10 22 g. With an interburst time At ~ 22 hr, 
we obtain M acc ~ 6.4 x 10 -9 M Q /yr. 

An extensive analysis of numerous bursts from 1636-536 was performed by 
Lewin et al. (1987), using EXOSAT data. Expansion was found in 3 bursts 
only, those ones with numbers 11, 26 and 27. Unfortunately no blackbody 
fitting during the expansion/contraction phase was performed for these bursts 
and the re-examination of archive data may provide valuable informations for 
this source. 

e) 1715-321 

Tenma observations (Tawara, Hayakawa & Kii 1984) have registered only 
one burst with expansion, and At is thus unknown. It may be that there 
were no bursts in the 3 months preceeding the burst under consideration, but 
it is also equally possible that At ~ 10 days, as it has been in 1979. The 
burst had T™ ax ~ 1.1 keV, T™{ n ~ 0.25 keV, corresponding to 7 = 0.8 and 
M env ~ 5 x 10 23 g. With At ~ 10 d, we obtain M acc ~ 9.2 x 10" 9 M Q /yr while 
for At ~ 3 months it is M acc ~ 10 -9 M Q /yr. The cooling time for such an 
envelope is t coo i ~ 10 days, so the situation may be fairly stationary for both 
values At considered. 

On July 20, 1982, Hakucho observed a very long (~ 300 s) burst with 
precursor, in which the luminosity stayed at its maximum value for almost 100 
s (Tawara et al. 1984a). The blackbody fitting gave T™ ax = 3.0 ± 0.2 keV, 
T™™ = 1.1 ± 0.1 keV. After this burst, in 29 days of observations, no events 
were detected (although the observations were not continuous, the effective 
observation time was ~ 120 hr). Generally saying, this source shows a weak 
activity, having produced only 3 bursts in 20 days in 1979 and a total of 4 
bursts in 155 days during 1979-82. Therefore, although the precise value of At 
for the burst we consider is unknown, a value ~ 10 d seems to be appropriate. 
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We obtain 7 = 2.29, M env » 2 x 10 23 g and the accretion rate before the burst 
was then M accr « 3.7 x 10~ 9 (10d/ At)M Q /yr, of the same order as with Tenma 
data. 

f) 1724-307/Terzan 2 
The only burst with expansion from this source was, probably, registered by 
Grindlay et al. (1980) during a 30 minute Einstein observation. The blackbody 
fitting of spectra produced T™^ x ~ 3.2 keV, T™\ n ~ 2 keV, resulting in the 
hardening factor 7 = 2.44 and M env ~ 5 x 10 22 g. The interburst time 
preceeding this burst is unknown, but after this event EXOSAT did not register 
any bursts in 12 hr of continuous monitoring. The value of accretion rate, 
M accr ~ 10~ 8 x (Id/ At)M & /yr may be equally sub- or super-Eddington, and 
no conclusion can be reached at the present stage. 

g) 1728-337 

The source was extensively observed by SAS-3 in 1976-78, and a total of 
60 bursts were registered (Hoffman et al. 1980). The average At was 3 -h 8 
hr. In only one of the bursts (on June 8, 1977) the expansion occurred. Before 
this event there were no bursts in a 56 hr period, although some bursts may 
have been missed because the Earth occulted the source roughly one-third of 
the time. The spectral fitting parameters were T™® x = 2.6 keV, T™™ = 1.7 
keV which give 7 = 1.98, M env ~ 4 x 10 22 g and, assuming the 56 hr value 
as the actual At, the accretion rate is M acc ~3x 10 _9 (56/ir/At) M Q /yr, just 
below the Eddington level. The thermal regime of the envelope between bursts 
is completely stationary since t coo \ ~ 1.2 hr. 

h) 1743-29 

The source is one of the Galactic Center X-ray bursters and it was observed 
with SAS-3. All the bursts reported by Lewin et al. (1976) had double and 
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triple peaks. The average At was ~ 1.5 d. Unfortunately, the spectral fitting 
was impossible due to the poor quality of the data. 

i) 1746-370/NGC 644 1 

Sztajno et al. (1987) detected 2 bursts with expansion, separated by 8.5 
hr, during a continuous 12 hr EXOSAT observation. The spectrum fitting 
parameters of both bursts are quite similar, T££ x = 2.2 keV, T™/ n = 1.2 keV. 
The application of our theoretical model yields 7 = 1.68, M env ~ 6.5 x 10 22 g, 
and an accretion rate M acc — 3.4 x 10~ 8 M Q /yr. The high accretion rate is in 
agreement with a cooling time scale ~ 10 4 s ~ At. 

The characteristics of the strongest bursts in 1746-370 seem to be not 
changing in time. In fact, a re-analysis of 2 bursts observed earlier with SAS-3 
(Li & Clark 1977) gives T™? x = 2.2, 2.4 keV, values quite similar to those 
derived from EXOSAT data. 

j) 1812-12 

The observations with Hakucho (Murakami et al. 1983) revealed two bursts 
with photospheric expansion, on 1982, August 18 and 22. The spectrum fitting 
produced T™ x ~ 2.5 keV; the data about the earlier phases of the burst are far 
poorer and do not allow the determination of T™\ n . Assuming, nevertheless, a 
value of T™™ ~ 1.2 keV as characteristic from other similar sources, we obtain 
7 = 1.9, M env ~ 9 x 10 22 g and M acc ~ 3.3 x lO" 9 (5d/At)M /yr. 

k) 1820-30/NGC 6624 
The source is located in NGC 6624, and this provides a reliable estimate for 
its distance of 6.4 ± 0.6 kpc (Vacca, Lewin & van Paradijs 1986). It exhibits a 
685-seconds periodicity, first discovered by Stella et al. (1987) with EXOSAT, 
and interpreted as an evidence for orbital motion. The orbital period is the 
shortest one known in LMXBs (see Parmar & White 1988 for a review) and is 
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consistent with a scenario in which the secondary is a low-mass, helium-rich 
degenerate star (Rappaport et al. 1987). Data on the time evolution of this 
period are somewhat controversial. It was observed to decrease over the years 
1976-1991 (Sansom et al. 1989; Tan et al. 1991; van der Klis et al. 1993a) with 
a time scale of ~ 10 7 yr, which might be caused by gravitational acceleration in 
the cluster potential or by a distant third companion in a hierarchical triple. The 
standard scenario, involving mass transfer from a Roche lobe-filling degenerate 
dwarf in an 11-minutes orbit around a neutron star, predicts a secular increase 
of the orbital period > 8.8 x 10" 8 yr _1 (Verbunt 1987; Rappaport et al. 1987) 
rather then a decrease. The most recent ROSAT observations (van der Klis 
et al. 1993b) do not provide any evidence of significant period decrease. It 
has been proposed that the secular variation observed in 1976-1991 could have 
been dominated by some changes in the position and shape of an occulting 
bulge on the disk rim. The phase shifts caused by this mechanism could, in 
principle, mimic the orbital period decrease, while the real period has been 
indeed increasing, according to the standard theory. 

As far as the evidence of photospheric expansion in this source is concerned, 
there are two series of observations which can provide information about the 
maximum and minimum color temperatures reached in the expansion phase, 
and the interburst time. The first is the series of bursts observed by Clark 
et al. (1977) with SAS-3 in May 1975 and March 1976 and the second is the 
sequence of 7 successive bursts with photospheric expansion observed by Haberl 
et al. (1987) with EXOSAT. The analysis of the latter data was presented in 
paper I and provided the unusually high values of M env ~ 9 x 10 23 g and 
M acc ~ 10- 6 M o /yr. 

Vacca, Lewin & van Paradijs (1986) reported 6 bursts with photospheric 
expansion in SAS-3 observations. The satellite had two independent detectors. 
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In some bursts both of them were providing enough information for the 
blackbody fitting, while, in other cases, only the data from one detector were 
of satisfactory quality. The estimates of the accretion rate are summarized in 
Table 1. The notation N/A stands for cases when T co i was not measured; the 
bursts are numbered as in Clark et al. (1977). In spite of being somewhat lower 
than the figure obtained from EXOSAT data, the present values M acc are always 
in excess of the Eddington value for He-rich matter, M Edd ~7x 10~ 9 M /yr. 
This confirms the suspicion that MXB 1820-30 is a very peculiar source which 
shows the highest mass transfer rate between X-ray bursters. Since bursters 
distinguish themselves among LMXBs because of high accretion rates, MXB 
1820-30 is, definitely, one of the observed sources with highest accretion rates. 
Together with the still unexplained abnormal behavior of the orbital period, 
our results make this source even more attractive for further investigations. 

I) 1850-087/NGC 6712 
The observations with SAS-3 revealed three bursts from the source, the 
first two separated by 17 hr and without signs of expansion, and the third, 
strongest one, with expansion (Hoffman et al. 1980). The precise value of At 
for the third burst is unknown, but a value ~ 1 d seems to be appropriate. 
The spectral fitting produced T™ ax ~ 2.7 keV, T™/ n ~ 1.6 keV. Then the 
hardening factor is 7 = 2.1, and M env = 5 x 10 22 g. The accretion rate is 
M accr ~ 9 x 10- g (ld/At)M Q /yr. Being close to the Eddington limit for He- 
rich matter, and only twice that one for a pure hydrogen envelope, such an 
accretion rate seems to be consistent with the fairly low persistent flux. 
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m) 1905+000 

Chevalier & Ilovaisky (1990) made a detailed spectral analysis of one burst 
with photospheric expansion detected during a continuous 19 hr EXOSAT 
observation. Although being non-perfect, the blackbody fitting gave T™ x = 
3.5 ± 0.7 keV, T™ n = 1.7 keV, resulting in 7 = 2.67 ± 0.53 and M acc ~ 
10~ 8 (ld/ At)M Q /yr. The thermal regime of the envelope should be stationary 
since t coo i ~ 10 4 s. 

n) 2127+119/M15 
One very strong burst with expansion was reported by Dotani et al. (1990) 
and Van Paradijs et al. (1990). Similarly to what was observed in 1608-522, a 
series of photospheric oscillations were detected with Ginga during the first ~ 30 
s of the burst in this source; this issue is addressed to in Lapidus et al. (1994). 
The blackbody spectral fit provided T c ™ ax ~ 2.8 keV, T™ n ~ 1.4 keV, although, 
like in some other sources, the spectra deviated from a blackbody in a way 
that cannot be described in terms of a spectral hardening alone (see discussion 
below). Applying our usual technique, we obtain 7 = 2.1 and M env ~8x 10 22 
g. The expectation time of a burst is unknown, but a value ~ few days seems to 
be appropriate. The estimated accretion rate M accr ~ 5 x 10~ 9 (3d/ At) M & /yr 
does not exceed, probably, the Eddington level. 
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IV. DISCUSSION AND CONCLUSIONS 

The results obtained in the previous section for 13 sources which exhibit 
photospheric expansion are summarized in Table 2. It can be seen from the 
table that, in several cases, the derived values of the accretion rate exceed the 
Eddington rate M Edd ~ 4 x 10~ 9 M & /yr for a "solar" composition. Such a result 
not necessarily contradicts standard steady accretion theories since, assuming 
spherical accretion, the total power released is ~ (GM*/.R*)M acc ~ M acc c 2 /6 
for a neutron star radius nearly three times the gravitational radius, as in our 
wind models. This means that accretion rates up to 6 times the critical one 
are still possible. This limit can be even higher if, as it seems more plausible, 
an accretion disk is formed around the neutron star. In disk accretion, in fact, 
only about half of the gravitational energy is released in the vicinity of the 
star surface (the rest being radiated away in the extended disk), and it is only 
this part of the energy release which can place a limit on the accretion rate. 
This argument would bring the upper limit for the permitted value of M acc to 
~ 10M Edd . Larger accretion rates are still possible in non-stationary accretion, 
for which the Eddington limit does not apply; in this case gravitational energy 
is temporarily stored in the accreted material in the form of internal energy. 

We note that, although our estimates of the accretion rate are model 
dependent, they do not imply any precise accretion scenario, that is to say 
they hold the same no matter how the gas is accreted onto the neutron star and 
where it comes from. It is apparent from table 2 also that the main source of 
uncertainty in the determination of M acc is At, which is often poorly known, 
while the evaluation of both 7 and M env appears to be more reliable. It seems 
reasonable to divide the sources in three classes, according to the estimated 
value of M acc : the range M acc ^ M Edd defines group I, M Edd ^ M acc ^> 10M Edd 
group II and M acc ^, ^0M Edd group III. Group I contains 4 sources: 1516-569, 
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1715-321, 1728-337 and 1812-12; in group II there are 8 sources: 0748-676, 
1608-522 (as observed with Tenma), 1636-536, 1724-307, 1746-370 (marginally 
belonging to group III), 1850-87, 1905+000, 2127+119. Out of the 13 sources 
listed in table 2, only 1820-30 definitely belongs to group III, confirming 
its peculiar nature, together with 1608-522 (as observed with Hakucho). In 
paper I it was suggested that the highly super-Eddington accretion rate in 
1820-30 could be reconciled with standard accretion scenarios on the basis 
of the unstationary nature of the accretion process. The comparison of the 
characteristic time needed to radiate away the gravitational energy release, 
which increases with M env , with the interburst time has shown, in fact, that, in 
this source, the envelope has no time to cool between two successive bursts, and 
therefore the persistent flux could not be directly related to M acc . The same 
argument applies to 1608-522 in the "high" state, as observed by Hakucho, when 
tcooi ~ 3At. It is remarkable that in all cases where M acc exceeds ~ lOMEdd, 
it is t coo i> At. This can be regarded as a confirmation of our model because 
the evaluation of the cooling time is completely independent on At. A further 
evidence is provided by the data of 1746-370, which seems to lie at the border of 
classes II and III and has M acc ~ lOMEdd, tcooi ~ At. Theoretical limits on the 
accretion rate for the onset of the helium thermonuclear runaway were placed by 
Fushiki & Lamb. According to their analysis of the stability of nuclear burnings 
in a static, solar composition envelope, a helium flash can be produced only if 
10~ 10 ^ M I (M Q /yr)^ 10~ 8 . The values we have derived are indeed within this 
range, with the exception of 1820-30 and, possibly, 1608-522, and this seems 
to be consistent with a picture in which unstable He burning is responsible for 
strong bursts with photospheric expansion. We note that, in the case of 1820- 
30, the accreted material is believed to be helium-rich, so that a straightforward 
application of Fushiki & Lamb's results to this source is not appropriate. As 



15 



far as 1608-522 is concerned, the accretion rate turns out to be well above 
10~ 8 M Q /yr only in one event, while other two observations provide values in 
the above range. The present estimates of accretion rates onto bursting neutron 
stars are also not in contradiction with the models of Taam et al. , who followed 
the time evolution of recurrent hydrogen thermonuclear flashes. They found, in 
fact, that accretion rates 10~ 10 M & /yr give rise to repeated H flashes while 
suggesting that unstable He burning may be produced only for larger M acc 
because the high transfer rate can inhibit convective mixing and increase the 
envelope cooling time. 

A definite improvement on the present method for estimating the accretion 
rate will come from self-consistent frequency-dependent transfer calculations 
which are now under way. They will provide, in fact, the hardening ratio for 
each wind model. This sequence of 7's will substitute the average constant 
value of 7 used so far, which was itself derived by matching spectral data with 
the models. The comparison of the computed hardening ratios with the present 
ones would also provide a further test for the consistency of our approach. 
We note, however, that the values of 7 listed in table 2 are quite reasonable. 
As we discussed in paper I about 1820-30, 7 ~ 2 is in good accordance with 
recent transfer calculations in expanding envelopes (Lapidus 1991) and, for all 
sources we examined, it is 1.5 < 7 < 2.7. The only exceptions are the Tenma 
observation of 1715-321, which give 7 = 0.8, while a "normal" value, ~ 2.3, 
was obtained from Hakucho's data. Frequency-dependent calculations will also 
allow the comparison of the self-consistently computed emerging spectrum with 
observations. This appears to be of particular relevance for strong bursts in 
which a pure blackbody fitting was often found to be unsatisfactory. Frequency- 
dependent radiative transfer calculations in expanding envelopes performed so 
far (Lapidus 1991; Titarchuk 1993; Turolla, Nobili, Zampieri & Lapidus, in 
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preparation), show, in fact, that a definite soft excess in addition to the general 
spectral hardening should be expected. 
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Table 1 



Accretion Rates in 4U/MXB 1820-30 (SAS-3 data) 



Burst # 


Detector 1 


Detector 2 




M acc (M Q /yr) 


M acc (M Q /yr) 


1 


1.7 x 10" 7 


9.7 x 10" 7 


4 


N/A 


1.7 x 10" 7 


6 


4.4 x 10" 7 


1.2 x 10" 6 


7 


2.8 x 10" 7 


N/A 


11 


N/A 


1.3 x 10" 7 


20 


1.4 x 10" 7 


4.2 x 10" 7 
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Table 2 



Accretion Rates and hardening ratios in X-ray bursting sources 



Source 


Mace 


7 


Notes 




(M Q /yr) 






0748-676 


<3x 10" 8 


1.7 


only lower limit for At available 


1516-569 


~10" 9 


1.7 


At only estimated 


1608-522 


8 x 10" 8 


2.7 


Hakucho, t coo i ^ At 




6 x 10" 9 


1.9 


Tenma burst G 




2 x 10 -8 


2.2 


Tenma burst J 


1636-536 


6 x 10" 9 


2.2 


Tenma burst H 


1715-321 


4 x 10 -9 


2.3 


Hakucho 




10" 9 -10- 8 


0.8 


Tenma, At unknown, drastic spectral softenin; 


1724-307 


~2x 10 -8 


2.4 


At only estimated 


1728-337 


<3x 10" 9 


2.0 


only lower limit for At available 


1746-370 


3 x 10 -8 


1.7 


t n r-. „ I r ^~ J ZAt 

"COOL 


1812-12 


~3x 10 -9 


1.9 


T m ) n unknown 

col 


1820-30 


10" 6 


1.5 


EXOSAT, He models, t cool ^>At 




10" 7 ^10- 6 


2.3 


SAS-3 (see table 1), He models 


1850-87 


9 x 10" 9 


2.1 


At only estimated 


1905+000 


~10" 8 


2.7 


At unknown 


2127+119 


~5x 10" 9 


2.1 


At unknown 
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